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Abstract 

Traffic models based on cellular automata have high computational 
efficiency because of their simplicity in describing unrealistic vehicular 
behavior and the versatility of cellular automata to be implemented on 
parallel processing. On the other hand, the other microscopic trafRc 
models such as car-following models are computationally more expen- 
sive, but they have more realistic driver behaviors and detailed vehicle 
characteristics. We propose a new class between these two categories, 
defining a traffic model based on continuous cellular automata where 
we combine the efficiency of cellular automata models with the accu- 
racy of the other microscopic models. More precisely, we introduce 
a stochastic cellular automata traffic model in which the space is not 
coarse-grain but continuous. The continuity also allows us to embed a 
multi-agent fuzzy system proposed to handle uncertainties in decision 
making on road traffic. Therefore, we simulate different driver behav- 
iors and study the effect of various compositions of vehicles within the 
traffic stream from the macroscopic point of view. The experimental 
results show that our model is able to reproduce the typical traffic 
flow phenomena showing a variety of effects due to the heterogeneity 
of traffic. 



1 



1 Introduction 



Traffic models are fundamental resources in the management of road net- 
work. A real progress in the study of traffic has obtained with the introduc- 
tion of the models based on cellular automata (CA j^. A cellular automaton is 
a collection of cells (sites) on a grid of specified shape (lattice) that evolves 
through a number of discrete time steps according to a set of local rules 
based on the states of neighboring cells, for more details see 0. CA models 
(CAMs) have the ability of being easily implemented for parallel computing 
because of their intrinsic synchronous behavior. A set of simple rules can 
be used to simulate a complex behavior, thus these models are conceptu- 
ally simple (for the theory and applications of CA as models for complex 
systems in several scientific fields, see [21 El [5]) • The traffic models based 
on CA are capable of capturing micro-level dynamics and relating these to 
macro-level traffic flow behavior. However, they are not as accurate as the 
usual microscopic traffic model^ such as the time-continuous car-following 
ones [1]. A basic one-dimensional CAM for highway traffic flow was first 
introduced by Wolfram, where he gave an extensive classification of CAMs 
as mathematical models for self-organizing dynamic systems [231 . Nagel 
and Schreckenberg proposed the first nontrivial traffic model (the NaSch 
model) based on CA for single-lane highway in 1992 [TTj. This is a time- 
and space-discrete model where the traffic road is divided into cells of 7.5 
m and it is defined on a one-dimensional array of a fixed number of sites 
with closed (periodic) boundary conditions. The velocity is expressed as the 
number of cells that a vehicle advances in one time step. The maximum 
velocity {vmax) is assumed to be 5 cells /sec (135 km/h) and every vehicle 
has the same target velocity Vmax- Each update of the movements of the ve- 
hicles is determined by four consecutive rules that are performed in parallel 
to each vehicle at each second. In literature, the CA traffic models are up- 
dated similarly where they only differ on the randomization step introduced 
to add some randomness, for a survey see [15| . Note that in all of them the 
cells represent the space as it is in the NaSch model, so in the sequel we call 
them as "NaSch-type" models [H HSl [201 [H] . 

In this paper, we propose a new traffic model via continuous cellular 
automata (CCA) which gets closer to the car-following models by introduc- 

^ Throughout this paper, abbreviation CA refers to both cellular automata (plural) and 
cellular automaton (singular). 

^We use the term "usual microscopic traffic models" to refer to the microscopic traffic 
models other than CA, which are in general defined by means of a system of differential 
equations. 
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ing some continuity without losing the computational advantages typical of 
CAMs. Continuous cellular automata (or coupled map lattices) are cellular 
automata where the states of the cells are real values in [0, 1], and the local 
transition rule is a real function. They have have been introduced by Kaneko 
as simple models with the features of spatiotemporal chaos, and have appli- 
cations in many different areas like fluid dynamics, biology, chemistry, etc. 
(for more details on couple map lattices, see [3110]). In this work we drop 
the restriction that the set of states is [0, 1], and we use the name CCA to 
consider CA where the set of state is in general an infinite set of the same 
cardinality of the continuum. Therefore, we consider a hybrid between the 
usual microscopic traffic models which are very accurate in predicting gen- 
eral traffic behavior but computationally expensive, and the CAMs which 
are very efficient due to their simplicity and intrinsic parallelism making 
them natural to be implemented for parallel computing. This process of 
passing from the typical coarse-granularity of CAMs to the continuity in 
space is done with a change of vision where we abandon the cell-space cor- 
respondence and assume that cells represent vehicles. Thus we consider a 
model with open boundary conditions where the number of cells is equal to 
the number of vehicles. In this way, we obtain the immediate advantage of 
having less cells to compute (just the number of vehicles). The continuity 
also gives us the possibility to refine the microscopic rules that govern the 
traffic dynamics, using fuzzy reasoning to mimic different real-world driver 
behaviors. All parameters of the decision process of the drivers are modeled 
individually by means of fuzzy subsets so that various types of drivers can be 
taken into consideration. Hence, we are able to study how the heterogene- 
ity influences the traffic macroscopically. The CCA model proposed here is 
defined first for a single-lane road and then we extend it to the multi-lane 
case. The assumption that cells represent vehicles makes this extension not 
as natural as it is in the NaSch-type models. 



2 The single-lane model 

In our CCA model, for the purpose of simulating different driver (vehicle) 
behaviors, we apply a multi-agent fuzzy system where the environmental 
information is transformed into the decision for the next time step action 
for each vehicle. Using this fuzzy logic-based system, we categorize the 
vehicles into types (kinds) where they have common characteristics such 
as the same perception of time or distance. In this way, we aim to study 
the effect of different composition of vehicles, i.e., the heterogeneity in road 
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traffic. Let us now describe our model for a single-lane road. Note that for 
the sake of simplicity, we fix the unit of time as one second. Consider the one 
dimensional CCA SC = {Z, 1^,^,5) where the lattice is the set of integers 
and the set of cell states S = (K xM+ xM+ xM x {L, 0, i?} x {L, 0, i?})U{_L}. 
A cell with the state _L represents a cell without a vehicle. The generic i-th 
non-empty cell is in the state ai{t) = {ki,Xi{t),Vi{t), Si{t),di{t),d'^{t)) where 

• ki represents the kind consisting of all the information (parameters) 
specified differently for each kind of vehicle, such as: the maximum 
velocity Vmaxi , the optimal velocity Vopti (the velocity which each kind 
of vehicle feels comfortable in traffic stream), the length li, the fuzzy 
membership functions, the maximum stress Smaxi-, the minimum stress 
Smirii, the probability functions of lane-changing to the right lane 
Pji.{x) and to the left lane Pi^{x). 

• Xi{t) is the position, defined as the distance from the origin of the 
road to the middle point of this vehicle, Vi{t) is the velocity, and Si{t) 
is the stress, a variable to keep track of how much the driver is above 
or below of his optimal velocity. The stress parameter is introduced 
to implement a more realistic driver behavior since we make the as- 
sumption that drivers usually tend to decelerate when they are moving 
with a velocity higher than their optimal velocity. We assume that if 
Si{t) < then the driver is stressed and wants to go faster, instead 
if Si{t) > then his stress derives from the desire to go slower. Note 
that the main usage of Sj(t) is related to the lane-changing process 
described in Section [3l 

• di{t) is the variable describing the desire for lane-changing to the left 
"L", to the right "i?" and staying on the same lane "0". d'^{t) is 
the variable showing from which lane the i-th vehicle is transferred: 
from the left lane "L", from the right lane "ii" or not transferred 
"0". Although these variables are used in the multi-lane model, the 
evaluations are done in the single-lane part. 

M IS a kind of one-dimensional extended Moore neighborhood defined by 
M{i) = {i — + l,i + 2), and (5 : S'^ — )■ S is the local transition function 
(local rule) acting upon a cell and its direct neighborhoods. The local rule 
is defined componentwise by 6{ai-i{t),ai{t),ai+i{t),ai+2{t)) = {ki,Xi{t + 
l),Vi{t + l), Si{t + l),di(t + l), 0). The variables Si{t + 1) and di{t + l) are de- 
scribed in Section[3l The space is updated as usual by Xi{t+1) = Xi{t)+Vi{t+ 
1) and the velocity by Vi{t + 1) = mm{vmaxi, ^xf {t),max{0, '^ii't) + 
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where Ai{t) is the acceleration calculated using the fuzzy decision modules 
and Axf{t) is the distance with the front vehicle from front bumper to rear 
bumper defined in ([T|). This front distance constraint in the choice of the 
updated velocity is introduced to make the model collision-free as in the 
NaSch model. However, since it is unrealistic, it is applied in our model 
only in the borderline situations where extreme decelerations are involved. 
Indeed, we have tried to avoid this constraint as much as possible by intro- 
ducing enough fuzzy rules to make the system more reactive to reduce these 
extreme situations. It is worth noting that in the case Ai(t) = 7.5 m/s'^, we 
essentially obtain the open boundary, space-continuous and deterministic 
version of the NaSch model |17j . 

Ai{t) is depending on the kind ki, the velocity Vi{t) and the variables 
defined as following. 



3. Next Front Distance (NFD): Ax+^(t) = Xi+2{t) - - ^ - |. 
This is introduced to have a better perception of the driver behavior 
of the next front vehicle since in reality drivers observe not only the 
vehicle just in front of them but also the vehicles ahead. 

4. Perceived Front Collision Time {PFCT): 



where Ci{t) is the parameter introduced for slowing down depending on 
the stress calculated by Qi{t) = ^'""^'(i)"'''*^ and T^{t) is Front Collision 
Time {FCT), the time that passes for the i-th vehicle to reach to 
(to collide with) the front vehicle, defined as T^(t) = ^.^^'^.^^^ ■ 



PFCT is a parameter which is a combination between FCT and an 
auxiliary time defined to keep the velocity close to the optimal velocity 
using a parameter depending on the stress variable. If the FCT is 
strictly negative, then the driver perceives an artificial vehicle with the 
FCT Ciit) to keep the vehicle close to its optimal velocity. Otherwise, 
PFCT takes the smallest FCT value between Cii't) and T^{t) to avoid 
collisions. 



1. Back Distance {BD): Ax~{t) = Xi{t) - Xi_i{t) - | - %i 

2. Front Distance {FD): 



Ax+(t)=x,+i(t)-x,(t)-^-|, 



(1) 





otherwise. 
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5. Worst Front Collision Time {WFCT): r+^(t) = ^|y^- The colli- 
sion time with the front vehicle in the case it suddenly stops. This 
parameter is introduced for safety reasons. 



6. Next Front Collision Time (NFCT): r+^(t) = ^.^t^Xj^^y The colli- 
sion time with the next front vehicle. It is clear that a vehicle never 
reaches to its next front vehicle, but NFCT is introduced to anticipate 
the braking maneuver of next front vehicle. 

7. Back CoUision Time (BCT): {t) = The collision time 
of the back vehicle. This is introduced to take into account the phe- 
nomenon where the back driver comes closer the front driver and forces 
him to accelerate. We call this situation as pushing effect. 

The approach of embedding fuzzy logic while dealing with a system de- 
scribed by continuous variables has been already introduced, indeed, there 
are several works based on fuzzy logic systems in car- following models [1]. 
In our model, we use a fuzzy rule-based system by means of fuzzy IF-THEN 
rules jl6j that determines the behavior of the vehicles in traffic stream. In 
this fuzzy system we introduce two fuzzy modules, formed by two sets of 
fuzzy rules, according to their significance. The first module has more im- 
portance since it is based on the information related to the front and the 
back vehicles while the second one has less importance since the rules in that 
module are related to the next front vehicle. As a consequence of having 
two different fuzzy modules, we have two outputs: Ai^i is the output coming 
out from the first fuzzy module and Ai^2 is the output coming out from the 
second fuzzy module. The first and the seconds sets of fuzzy rules used in 
the simulator are summarized in Table [1] and Table [21 respectively. After re- 
ceiving the inputs r+p(t), r+^(t), -r^^rCt), T~{f), Axf{t), Ax+^(t), Axfit), 
Vi{t) determined by the environment, the final decision of the acceleration 
Ai(t) is returned by the function (to avoid cumbersome notation, we omit 
the dependency of t): 



FiAi^uA 



i,2) 



mm{Ai^i,Ai^2) if < 0, 

^^^4^ if^,,i >0AAi,2 < -0.25, (2) 

Ai 1 otherwise. 



where we give more weight to the decision taken by the first module. Indeed, 
in the simulation, we have noticed that without these kind of constraints 
the vehicles were tending to slow down too much. This is also the reason of 
splitting the fuzzy module into two parts. Moreover, the second module and 
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consequently the second acceleration gains significance only in the case of 
emergencies such as a sudden breakdown or deceleration of the next front ve- 
hicle. The linguistic terms (properties) used for the acceleration outputs Ai 
are: PB=Positive Big, PM=Positive Medium, PS=Positive Small, Z=Zero, 
NS=Negative Small, NM=Negative Medium, and NB=Negative Big. 



Table 1: First fuzzy set of rules 





R'ont Distance (FD) 


Velocity (vj) 




BIG 


MEDIUM 


SMALL 


VERY SMALL 


SMALL 




BIG 


(i.j NOT SMALL) 


(ii. NOT SMALL) 


ZERO 


ZERO 


PB 


PFCT 


PM 


PS 






(Jam Situation) 




MEDIUM 


ZERO 


ZERO 


NS 


NS 






SMALL 


NM 


NM 


NM 


NM 






VERY SMALL 


NB 


NB 


NB 


NB 







WCT 1 


SMALL 1 


NS 


1 NM 1 


NM II 



In the first fuzzy set of rules we have also the ones related to the back 
vehicle, thus we consider the pushing effect where BCT is "Very Small" 
and BD is "Very Small". In this case, we apply the following rule formed 
by taking into consideration the safety with the front vehicle: (IF PFCT is 
BIG AND FD is BIG) OR (IF PFCT is BIG AND FD is MEDIUM) OR 
(IF PFCT is MEDIUM AND FD is BIG) OR (IF PFCT is MEDIUM 
AND FD is MEDIUM) then the acceleration choice is PS. 



Table 2: Second fuzzy set of rules 





Next Front Distance (NFD) 


VERY SMALL 


SMALL 


MEDIUM 


BIG 


NFCT 


VERY SMALL 


NB 


NB 


NM 


NM 


SMALL 


NM 


NM 


NS 


NS 


MEDIUM 


NS 




BIG 


NS 



These rules are based on some common sense of driver behaviors in- 
cluding some experiences. Although the rules are the same for each kind 
of vehicle, they have different "weights" depending on the definition of the 
membership functions of different kinds. In this way, it is possible to give a 
description of a variety of behaviors, such as the behavior of a long vehicle 
driver or a driver with low reflexes. For instance, a person with low reflexes 
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perceives a time of collision of 5 seconds as a very short time, however a 
person with higher reflexes probably would feel comfortable with that time 
of collision. 

After the fuzzy system receives all the inputs, we first determine the 
degree of membership of each input. Then by taking the minimum value of 
the images of the inputs (corresponding to the fuzzy logic AND operation) 
in each rule, we obtain the weights of the fuzzy rules. This evaluation is 
described in general terms as following. 

Let BI, ■■■,Bj^_ and be fuzzy subsets with the membership functions 

fi^j , . . . , fi^j and fi^j , respectively, and let be the fuzzy rules defined 
as: ^ 

: IF x{ is b{ AND . . . AND is B^^ THEN y is & 

with 

lij^j {x{) A ... ^^iJJJ (x{.) = min <^ ^ , . . . , /i^, 

for 1 < j < m, where m is the number of fuzzy rules. 

In our fuzzy system, there are also rules including the form: xi is not Br 
with the degree of membership fJ.^^j {xi) for some r € These rules 

which are the first two rules of the first fuzzy decision module, are formed 
with the purpose of emphasizing the vehicle is not in a jam situation (see 
Table [H where Vi is not small). For an input xi, the degree of membership 
of not having a property Br ("not being i?r") can be written as: 

The output of the previous process so far is a fuzzy set, so we should con- 
vert our fuzzy output set into one single number as the output of the fuzzy 
system, which is the "acceleration" in our case. Recall that since we have 
two fuzzy modules, there are two outputs of the system: Ai^i and Ai^2 which 
are used to evaluate Ai{t) at time t by the function F, see Equation [2j The 
conversion process is called as defuzzification and there are many defuzzifica- 
tion techniques to obtain a crisp value, such as the center-of-gravity (COG) 
and the weighted average formula (WAF) method. In our model, since we 
would like to use a simple method that does not require too much compu- 
tational power, we define a new method and call it as generalized weighted 
average formula (GWAF), i.e., we generalize the WAF to the case where the 
membership functions are not necessarily symmetric. We describe it now in 
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general terms. Suppose that each j-th rule receives the values x{, ■■■,xl,, as 
inputs. Let 

= minj^^, {x{), ... (xyj 

be the weights for each rule j and let = fi~](w^) be the preimage of the 
weight of j-th rule. We assume that fj.Qj does not have any plateau, since 
we do not want to contain intervals (this is done to have a discrete set). 
The defuzzified output is thus calculated by: 

3 The multi-lane model 

In this section, we extend the single-lane model to the multi-lane case. This 
extension is not trivial as it is in the NaSch-type models where adding a 
lane simply means adding an array of cells and where the local transition 
function can naturally be extended. This is a consequence of having a clear 
physical interpretation of the model given by the fact that space is rep- 
resented by cells. However, as a consequence of not having the cell-space 
correspondence, first of all we do not have the natural order which makes the 
extension of the local transition function so easy to achieve, and secondly 
not all the configurations c : Z ^ E of our model represent a physical situa- 
tion. For this reason, we consider just the configurations that have physical 
meaning, i.e., two adjacent cells n, n + 1 are in the states where the positions 
fulfill the constraint Xn{t) < We denote the set consisting of such 

configurations by Confp{SC). The support of a configuration c G Confp{SC) 
is a maximal interval of integers with i < j such that for all A; < i and 
k > j, c{k) =±. 

The union of two arrays of cells representing a road with two lanes is 
a natural candidate for a multi-lane model. Therefore, we first present 
our multi-lane model as a union of interacting single-lane CCA where the 
interaction is given by a transfer operation, and then we prove in Prop. 13.11 
that this model can actually be simulated by a CCA. 

The process of lane-changing behavior consists of several steps. First the 
driver feels the desire of lane-changing, then if the conditions of the place 
in the lane where the driver wants to move respect some safety constraints, 
the driver can finally perform the maneuver. Regarding the decision of 
performing a lane-changing, we first check if the i-th vehicle with the state 
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ai{t) = {ki,Xi{t),Vi{t),Si{t),di{t),d^{t)) desires to change lanes, i.e., di(t) ^ 
0. The desire of lane-changing is depending on the stress variable Si{t). 
Having positive stress means that the driver is above of his optimal velocity, 
thus he tends to change lane to the right (slower lane). Otherwise, the driver 
is nervous and he desires to go faster, thus he tends to change lane to the 
left (faster lane). The stress cannot be increased (or decreased) arbitrarily, 
so we define the two parameters Smim and Smaxi which are the maximum 
amounts of negative and positive stress that a driver can tolerate. Thus 
if the updated stress Si{t + 1) exceeds the boundaries Smm, and Smaxi 
simply assign it to these limit values. The stress is updated by 

r Sacc{t)/2 if r+(t) < and ^ < Sacdt) < 0,(3) 

S^(.t + 1) = I Saccit) •(! + $) if T+ (t) > and ^ < Saccit) < 0,(4) 

I Sacc{t) otherwise. 

where the accumulated stress is defined by Saccit) = Si{t) + {vi{t + l) — Vopti)- 
X(t) and X(t) is a random variable distributed uniformly, X(t) ~ U{0, 1). 
Condition ([3]) is used to avoid frequent lane-changing in the case of a jam 
situation, if the queue is moving. We also embed the strategy of trying to 
change lanes instead of braking in the case the front vehicle is close and 
tends to brake. We model this effect with the condition (jj]) considering a 
factor <I> representing how much the i-th. vehicle is in the situation where 
the front vehicle is close and going much more slower than the i-th vehicle. 
$ is calculated by the membership value of (Tj^(i) is VERY SMALL A 
Axf{t) is MEDIUM)v{T+{t) is VERY SMALLAAxf{t) is SMALL)V 
(r+(t) is SMALLAAxf{t) is MEDIUM)^ {t+ (t) is SMALLAAxf{t) is 
SMALL). It is used to increase the amount of stress and consequently to 
have more probability of changing lanes. 

We now describe the stochastic process that evaluates the desire of lane- 
changing for the next time step di{t+l) = Eval(^jr -ji^{ki,Vi{t), Si{t)) described 
in Algorithm [H The decision of this action is made by means of a Bernoulli 
process B{2,p). The probabilities of changing lane to the left and to the 
right, which are in general different, are calculated by two functions Pl.{x) : 
[0, 1] [0, 1] and P{i.{x) : [0, 1] — ?• [0, 1] contained in the kind ki. The reason 
is that some kinds of vehicles tend to move more to the left lane while some 
others tend to move more to the right. For instance, long vehicles prefer to 
go to the right lane more than left, i.e., they tend to stay more on the right 
lane. The variable used to calculate such probabilities is the normalized 
stress nsi{t) = m.ax.{si{t) / Smim, Si{i) / Smaxi)- Thus if Si{t) > then the 
choice to stay or go to the right is done by B{2, Pji.{nsi{t))), otherwise 
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Algorithm 1 The pseudo-code for evaluating Eval 



procedure Evalf^ jif{ki,Vi{t), Si{t)) 
if Si{t) > then 

Execute the Bernoulh Trial X ~ B(2, Pr. {nsiit)) 
if X = 1 then 

di{t+l) = R 
else 

di{t+l) = 
end if 
else 

nsi{t) = Si{t)/s 

mini 

Execute the Bernoulli Trial Y ~ B(2, PL\nsi{t)) 
if Y = 1 then 

Execute the Bernoulli Trial Z ^ B{2, ij.veis{vi{t))) 
if Z^l then 
if 7^ = then 

d,(t + l)^L 
end if 

if £ = then 

d,{t + l) = R 
end if 

Execute the Bernoulli Trial W ~ S(2, 0.7) 
if W = 1 then 

d,(t + l) = L 
else 

d^it + l) = R 
end if 
else 

d,{t+l) =L 
end if 
else 



end 
end if 
end procedure 



diit + 1) = 



by B{2,PL^{nsi{t))) to go to the left. If the latter process is resulted as 
going to the left, we take into consideration also the jam situation where the 
drivers randomly move to the left or right to find an emptier lane. For this 
purpose we use the velocity parameter and we need the fuzzy set of "Velocity 
SMALL" with the membership function denoted by Hveis to perform this 
evaluation which is done by the Bernoulli process B{2, fiyeis{vi(t))). If the 
result of the outcome is positive, meaning that the traffic is indeed jammed, 
the decision for lane-changing depends on the relative position of the lane, 
i.e., if the vehicle is in the most left (right) lane then it can only move to 
the right (left). For this reason, the single-lane model SC^^cn) depends on 
two parameters C,TZ which are Boolean variables used to describe whether 
or not there exists a lane on the left or on the right, respectively. Thus if 
there is no lane on the left (right) of the driver, i.e., C = Q {IZ = 0), he can 
only move to the right (left) lane. If there are lanes both on the left and 
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right, the choice is obtained using the Bernouhi process fi(2,0.7) where the 
decision is "moving to the left lane" with the probability of 0.7 and "moving 
to the right lane" otherwise. These probabilities are decided to be different 
since we make the assumption that the drivers desiring to go faster usually 
tend to move to the left lane more than the right even in a situation of traffic 
congestion. 

After the desire of lane-changing, we perform the transfer of a vehicle. 
The transfer clearly changes the configuration of the single-lane CCA, thus 
we need to introduce some operations to describe this process. Given a vehi- 
cle represented by a state a = (/c, x, v, s, d, d') € S \ {-L} and a configuration 
c € Confp{SC), the inserting operator at position n G Z is the function 
Insn '■ S \ {_L} X Confp{SC) — > Confp{SC) that changes the configuration 
from c into Insn{cr,c) by shifting the states of all cells one step starting 
from the n-th position and setting the state of the n-th cell to the value a. 
The right-inverse of InSn is the deleting operator at position n denoted by 
Deln : Confp{SC) — )■ Confp{SC). The index n gives the position where we 
insert or delete the content of a cell. The function Insn{(y,c) is well defined 
only if the insertion of the vehicle still generates a physical configuration. 
However, before the insertion, we check what is the relative position of the 
vehicle that wants to perform a lane-changing with respect to the lane that 
it is going to enter. For this reason, we need to define an index operator 
Indx : S \ {_L} x Confp{SC) — )> Z such that Indx{a, c) represents the index 
where the vehicle with state a would go if we would try to insert it into the 
configuration c. 

In our model, the process of lane-changing is based on some safety crite- 
rions which check the possibility of executing a lane-changing by considering 
the situation in the target lane. These criterions, obtained from some simu- 
lation experiments, guarantee that after the lane-changing there will be no 
danger (avoidance of collision) with the back and the following vehicle on 
the target lane. Consider the vehicle represented by a that wants to enter 
in a lane with configuration c, thus the corresponding position in the target 
lane is between the j — 1-th and j-th vehicle with velocities vj-i and Vj 
where j = Indx{a, c). Suppose that the back and front distances of this 
entering vehicle with respect to the j — 1-th and j-th vehicle are denoted by 
A~ and A"*", respectively, then the safety constraints are: 

A" > max(0, vf3i - v + \vj-i - v\ + 3), A+ > max(0, v^'^^ - vj + 3) (5) 

The condition for performing a lane-changing to the left (right) lane can 
be summarized by the Boolean variable transL{a,c) {transR{a,c)) which 
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consists of the safety criterion ([5]), the desire to change lane to the left 
(right), i.e., d = L {d = R) and not being already transferred, i.e., d' = (to 
avoid multiple transfers in one step of simulation). In the process of lane- 
changing, we also need to keep trace where the vehicle comes from because 
the transfer of a vehicle is seen as a process consisting of copying and erasing 
steps. For this purpose we define the copy of a as: 



a 



cp 



{k, X, V, s/5, L, R) if cr = {k, x, v, s, L, 0), 
{k, X, V, s/5, R, L) if cr = {k, x, v, s, R, 0). 



where we decrease the stress from s to s/5, firstly because we make the 
assumption that when a vehicle changes lanes there is a sense of satisfaction 
reducing the stress, and secondly this change helps to reduce the ping-pong 
phenomenon [12^ [20] of vehicles continuing to move from one lane to another 
in case of traffic congestion. If a vehicle with state a desires to go to the left 
(right) lane, then we copy a into c, i.e., we update c into a new configuration 
denoted by a c (a c)€ Confp{SC) defined by the following equation 

{Xe{L,R}y. 



Insj{a'^P, c) if trans x{cr, c), where j = Indx{a, c) 
c otherwise. 



Besides, we need an operator to update a configuration c' in which we have 
to erase the state that has already been copied into c. The erasing procedure 
is accomplished by changing c' into the new configuration c'\uj G Confp{SC) 
defined as: 



c'\uj 



if (c'(i) = {h, y, w, 5r, R,0) A u = (h, y, w, r, R, L)) 
Deli{c') \/{c'{i) = {h, y, w, 5r, L, 0) A a; = (/i, y, w, r, L, R)), 

where i = Indx{u}, c') 

d otherwise. 



We extend these functions to the binary operations and \ on 

Confp{SC) to describe the transfer process of all the vehicles from one 
lane to another. We call and as the left and right copying op- 
eration, respectively, and \ as the erasing operation. The extension is made 
in the following way. Let c, c' G Confp[SC) be two configurations with 
c(z) = oJi /_L for i = m, . . . , M and c'{i) = ai 7^_L for i = n, . . . , A^, where 
[m, M], [n^N] are the supports of c, c', respectively. We define the config- 
urations c' c inductively as follows: let eo = c and e^ '■= CTk '^L ^k-i 
for = 1, . . . , — n + 1 then we define c' c := e^-n+i (analogously 
for ^r). On the other hand, for c' \ c, let = d and ■= gk-i \ fo^ 
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Algorithm 2 The pseudo-code for one time step evolution of the multi-lane 

model 

1: procedure Update(co, . . . , ca/-i) 



2: for i = M ~ 1 do 

3: if i = then 

4: Ci ~ (co ci) 

5: Co ~ Co \ (co ci) 

6: end if 

7: if < i < M - 1 then 

8: Ci_i — (ci Ci_i) 

9: Ci ■— d \ (ci )-^L Ci_i) 

10: Ci+i (ci Ci+i) 

11: Ci ■- Ci \ (ci Ci+i) 

12: if i = 1 then 

13: Co ■— 5(*oa)(co) 

14: else 

15: Ci_i := 5('i,i)(ci-i) 

16: end if 

17: end if 

18: if i = A/ - 1 then 

19: CA/-2 ~ (ca/_i CA/-2) 

20: cm-1 := ca/_i \ (ca/-i CM-2) 

21: CAf_2 := (5{*i i)(ca/-2) 

22: CAf_i := q)(ca/-i) 

23: end if 
24: end for 
25: end procedure 



k = 1,...,M — m + 1 then c'\c:= gM~m.+i- Using these operators it is easy 



to see that the process of the transfer of vehicles from c' to c, where c is the 
configuration of a lane on the left of the lane with configuration c', can be 
easily described by transforming c into c' c and c' into c' \{c' c). 

We now present our multi-lane model using these copying and erasing 
operations. Suppose that there are M > 2 lanes on a highway. We model 
each lane using the single-lane CCA model SC. Recall that SC depends on 
the parameters C,TZ, thus we rewrite it as SC^^cn) = (^i ^j-^j ^(£,7^))• 
this way, we can associate to the road the ordered M-tuple: 

where if M > 3 we have M — 2 copies of 5>C(i i), and where 5£(o,i) and 
5/^(1 0) represent the left- and the right-most lane, respectively. In the case 
M = 2, we consider just the pair (5>C(o,i)5 '^^(1,0))- Suppose that these M 
number of CCA are in the configurations (cq, . . . , cm-i) G Confp{SCY'^ . In 
our multi-lane model, we scan each lane and we transfer the vehicles to the 
adjacent lanes. After this process, for each lane we apply the single-lane 
CCA model to update the configuration and this update is done by means 
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of the global transition function of SC(^a,b) denoted by (^(^q^), a,b & {0)1}- 
In this way, we have a new array of configurations (cg, . . . , c'pj_-^), and this 
process represents 1 sec. of simulation. Moreover, the order with which the 
transfer is performed is from the left-most lane to the right-most one, and 
this is done to satisfy the precedence requirement in European roads. In 
Algorithm[2l it is described the updating process Update : (cq, . . . , cm-i) i-^ 
(cq, . . . , c^_i). We now show that Algorithm [2] can be simulated by a CCA 
which implies that our multi-lane model is framed as a CCA model. 

Proposition 3.1. There is a stochastic CCA multi-lane model which sim- 
ulates Algorithmic 

Proof. We define the stochastic CCA: 

MC = {Z,n,M,A) 

where 

• = {Confp{SC) X {copy , erase] x N^) U {-L}, where _L is the state 
associated to the empty cells representing cells with no lane. 

• M{i) = (i — 1, i, i -|- 1) is the von Neumann neighborhood. 

• A : — 7> O is the local transition function such that 

A(cj_i,wo,t^i) = 
defined in Algorithm [3] with inputs 

ujj = {cj,Xj,Mj,Pj,Kj), j = -1, 0, 1. 

If we consider M lanes with the configurations cq, . . . , cm-i, we associate to 
M£ = {Z,Q,M,A) the configuration 

C = {U}Q, . . . ,U}M-l) 

where = (cj, copy, M, i, 0) for i = 0, . . . , M — 1. It is easy to see that 
applying 2M times the global transition function A* to C, we obtain a new 
configuration 

A*^^^(C) = (a;^,...,^,_i) 
with = {c[, copy, M, i,0) for i = 0, . . . , M — 1 where 

(c'o, . . . , 4/„i) = Update{cQ, cm-i), 
and Update{cQ, . . . ,cm-i) is the function defined in Algorithm [2J □ 
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Algorithm 3 The pseudo-code to compute the local transition function A. 



1: 
2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 
10 
11 
12 
13 
14: 
15: 
16: 
17: 
18: 
19: 
20: 
21: 
22: 
23: 
24: 
25: 
26: 
27: 
28: 
29: 
30: 
31: 
32: 
33: 
34: 
35: 
36: 
37: 
38: 
39: 
40: 
41: 
42: 
43: 
44: 
45: 
46: 
47: 
48: 
49: 
50: 
51: 
52: 
53: 
54: 
55: 
56: 
57: 
58: 
59: 
60: 
61: 
62: 
63: 



procedure A (ixJ-i, uju, lji) 
if LJo —1. then 

uj'o =± 
else 

it Ko = then 

if Xo — copy then 

if Po = Kq + 1 then 

Co ~ (C-l Co) 

end if 

Xo := erase 
else 

if Po = Ko then 

Co := (co \ ci) 
end if 

Ko := A'o + 1 mod Mq 
Xo ~ copy 
exit 
end if 
end if 

if < Ko < Mo - 1 then 
if Xo = copy then 

if Po = Ko ~ 1 then 
Co := (ci >-*L Co) 
end if 

if Po = Ko + 1 then 
Co := (c_i >->fl Co) 
end if 
Xo := erase 
else 

if Po = Ko then 

Co := (co \ c_i) 
Co := (co \ ci) 
end if 

Ko := Ko + 1 mod Mq 
Xo ~ copy 
exit 
end if 
end if 

if Ko = Mo - 1 then 
if Xo = copy then 

if Po = Ko - 1 then 

Co := (ci >-^L Co) 
end if 
Xo :— erase 
else 

if Po = Ko then 

Co := (co \ c_i) 
end if 

if Po = then 

Co := Slo,i){co) 
end if 

if < Po < Mo - 1 then 

CO := '5(i_i)(co) 
end if 

if Po = Mo ~ 1 then 

CO := '5(i_o)(co) 
end if 

Ko := Ko + 1 mod Mo 
Xo ~ copy 
end if 

end if 16 
end if 
end procedure 



4 The experiments and the simulation results 



The model presented here is a CAM which is intrinsically parallel. We first 
implement the model using the programming language Python since it is 
a high level language making the implementation of the model faster and 
easier. Indeed, during the phase of the development of writing the code, we 
use an object-oriented philosophy, especially while passing from the single- 
lane to multi-lane case. Thus, it is easier to modify and to rewrite some 
parts of the code to tune it better. Using this code we run a series of 
experiments for a first evaluation of the model. Part of the data obtained 
by these experiments have been presented in the EWGT 2012 conference 
[26], and the complete version is contained in the PhD thesis [25]. We also 
implement the model using CUDA to take advantage of the power of the 
modern graphic cards to perform parallel computation. Indeed, the single- 
lane model can be naturally parallelize simply by giving a cell (a vehicle in 
our case) to each thread of the GPU. In Section [5] we give a brief description 
of this process. 

4.1 Setting the kinds of vehicles and the experiment scenar- 
ios 

In this section, we describe the setting up conditions for the kinds of vehicles 
used in a series of simulations. In our experiments, we use just two kinds of 
vehicles which we call as passenger vehicles and long vehicles. The vehicles 
have the following parameters: 

• passenger vehicles: Vmax = 36 m/s, Vopt = 28 m/s, the length 4 m, 
the natural acceleration noise: 0.2 m/s2 (see [7]), 

Smax — 500 m, 

Smin = —450 m, the function of the probability of lane-changing to the 
right lane Pr{x) = x, the function of the probability of lane-changing 
to the left lane Pl{x) = x. 

• long vehicles: Vmax = 25 m/s, Vopt = 20 m/s, the length 9 m, the 
natural acceleration noise: 0.1 m/s^ (see [^), Smax = 300 m, Smin = 
—700 m, the function of the probability of lane-changing to the right 
lane Pr{x) = x, the function of the probability of lane-changing to the 
left lane Pl{x) = x^-'^^. 
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Figure 1: The membership functions of the passenger vehicles for the front 
distance 

Regarding the fuzzy membership functions, we tune them according to a 
questionary posed to a group of drivers. The results represent how much the 
claims are related to each interviewee's perception while driving. We then 
consider the corresponding membership functions by interpolating linearly 
these data. However, in the simulations, we do not use these functions 
based on the data. Instead, we approximate them with the triangular and 
trapezoidal membership functions which are more common and easier from 
the computational point of view. We notice from the experiments that this 
kind of approximation do not alter the results of the simulation significantly. 
Firstly, because the sample is too small so the data obtained is subject to 
errors, and secondly, the obtained membership functions are not so different 
from their approximated versions and the simulation can gain a boost of 
a factor 1.5. In Figure [U we see one of the membership functions of the 
passenger vehicles for the perception of front distance adjusted according to 
the collected data and the approximated version, respectively. 

For the first evaluation of our model, we simulate a piece of highway 
near an off-toll plaza. The parameters used in the simulations are: 

• Length of the road L 

• Number of lane M 

• Number of iterations (the simulation time) 

• Number of repetitions of the same experiment 

• Emission rate A: the number of vehicles entering this piece of highway 

• Influence radius p of the off-toll plaza, note that p = —1 means that 
there is an open road tolling 
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• Obstacle: An obstacle can be placed on the left-most lane or on the 
right-most lane 

We model the entrance of the each lane using a Poisson stochastic process. 
The probability of emitting at least one vehicle in a lane in the time interval 
[0,t] is calculated by a Poisson distribution. If the random test succeeds, 
i.e., if the result is "emit a vehicle" then we randomly choose a kind of 
vehicle where the choice is a "passenger vehicle" with the probability of 
1 — p and a "long vehicle" with the probability of p, and in the experiments 
we consider p = 0%, 10%, 20%, 30%. The emission rate is set to X/M for 
each lane where we consider A = 0.25,0.5,1,1.5,2 veh/s. For instance, if 
A = 1.5 and M = 3, the average number of vehicles entering the piece 
of highway is 90 veh/min which means that each lane at the beginning of 
this road is charged of around 30 veh/min corresponding to a situation of 
heavy traffic since the maximum capacity of a lane is considered around 
40 veh/min [13]. Regarding the off-toll plaza, we introduce a parameter 
called influence radius which is the interval of the road segment in which 
the vehicles are captured when they are exiting. This parameter is used to 
control the processing time of each vehicle which influences the throughput 
(the number of vehicles processed in every 10 seconds). The influence radii 
considered are p = 10, 25, 50 m and "-1" which means that the off-toll plaza 
is not visible to the vehicles, thus they simply act as if there is no any off- 
toll plaza. This situation corresponds to the case of an open road tolling 
payment system where the vehicles do not need to slow down and stop to 
make payment. Every 10 seconds the simulator checks the number of vehicles 
processed (exited) and it calculates the average latency, i.e., the average of 
the time spent to travel from the entrance to the exit. We also place an 
obstacle in some experiments either on the left- or on the right-most lane 
with a dimension of 2L/5 in the middle of the piece of highway to analyze 
a bottleneck phenomenon. In the simulations, each experiment has 1000 
iterations and it is repeated 100 times. The execution time takes around 
6.6 hours using one core of a computer equipped with a 16 core Xeon at 2.7 
GHz (X5550) with 16 GB of RAM running Debian Linux (kernel 2.6.32). 

4.2 Analysis of the experimental results 

The macroscopic parameters considered for the purpose of analyzing the 
general behavior of traffic are evaluated as following. The density describing 
the number of vehicles per unit length of the piece of highway (measured 
in vehicles per meter) at time t is calculated by D(t) = where N(t) 
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is the total number of vehicles at time t and L is the length of the road 
representing a piece of highway (5 km long) near the off-toll plaza. The 
number of vehicles is not constant since the entering rate of the vehicles 
(depending on a Poisson stochastic process) is different from the exiting 
one, so we obtain different densities. In our experiments, the density can 
reach at 0.25 veh/m/lane as maximum since the length of a passenger vehicle 
is assumed to be 4 m long. The average velocity, i.e., the averaged sum of 
the velocities at time t is calculated by Vav{t) = "^i^Ji Vi/N{t), and finally, 
the flow is evaluated by q{t) = D{t) ■ Vavit)- 

The analysis of traffic flow is typically performed by constructing the 
fundamental diagram (the flow-density diagram) that determines the traffic 
state of a roadway by showing the relation between flow and density. There- 
fore, the computations described above are used to plot the fundamental di- 
agrams where we analyze different traffic phenomena, see for instance Figure 

El 

In the simulation results, we notice that the different composition of 
vehicles obtained by changing the percentages of the long vehicles is an im- 
portant factor in influencing the flow. It is seen in Figure [2] that adding even 
a small percentage of long vehicles to the traffic stream changes the funda- 
mental diagram significantly. The heterogeneity effects also the throughput 
as it is predictable since long vehicles are slower and so in the queue near 
the off-toll plaza it takes more time to move and to get processed. 

The experiments show that the model is able to reproduce the typical 
traffic flow physical phenomena such as the three phases of traffic flow |llj : 
Free flow, synchronized flow and wide-moving jam, see Figure [3j Free flow 
corresponds to the region of low to medium density and weak interaction 
between vehicles. In general, the slope of the fundamental diagram in the 
free flow phase is related to the speed limit, meaning that in this phase 
vehicles can move almost at the speed limit. Instead, in the free flow phase of 
the fundamental diagram in our experimental results, the slope is related to 
the optimal velocity, meaning that in this phase the vehicles can move almost 
at the optimal velocity. The reason is that in our model the vehicles do not 
aim to reach to the maximum velocity as it is in the NaSch-type models, 
but they tend to go with their optimal velocity. Free flow is characterized 
with a strong correlation and quasi-linear relation between the local flow 
and the local density [19]. The synchronized flow presents medium and 
high density while the flow can behave free or jammed. In other words, it 
is defined by the interaction between the vehicles and is characterized by 
an uncorrelated flow-density diagram. However, this phase is not clearly 
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Figure 2: Comparison of fundamental diagrams with 3 and 4 lanes, respec- 
tively 
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Length of the road: 5000m, Number of Lanes: 3, Emission Rate: 2 vch/sec, Inf Radius: 10 
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Figure 3: Traffic phases in the fundamental diagram: Free flow, synchronized 
flow and wide-moving jam, and the cross-covariance between the flow and 
density 



understood in the context of CA and not observed in most of the NaSch- 
type models. It probably requires the presence of sources (on-ramps, on-toll 
plaza) and sinks (off-ramps, off-toll plaza), see |22j . The wide-moving jam 
phase represents the situation where the traffic is jammed (congested). In 
this phase, an increase in the density results with a decrease in the flow. Let 
us now consider the cross-covariance between the flow q(t) and the density 
D{t), see Figure [3l 

. j^. {q{t)D{t))-{q{t)){D{t)) 

where the brackets (•) indicate averaging the values obtained in all the ex- 
periments at time t. In the free flow phase, the flow is strongly related to 
the density indicating that the average velocity is nearly constant. For large 
densities, in the wide-moving jam phase, the flow is mainly controlled by 
density fluctuations. There is a transition between these two phases where 
the fundamental diagram shows a plateau when the cross-variance is close 
to zero. This situation with cc{q, D) ~ is identified as the synchronized 

flow [laiiH]. 
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It is seen in Figure [3] that we reproduce these relations where the free 
flow phase initiates with the situation of cross-covariance close to 1 and 
continues with a positive cross-covariance, the synchronized flow phase is in 
the region where the cross-covariance is close to zero and the wide-moving 
jam phase corresponds to the situation where the cross-covariance is negative 
(anticorrelation) . 

The plateau formation has the dependency also on the throughput. In 
Figure HI it is seen this situation where we have compared the fundamen- 
tal diagrams for different influence radii which are related to the average 
throughputs. Observing any of these three plots, we see that when the aver- 
age throughput is increased, there occurs a more immediate transition from 
the synchronized flow phase to the wide-moving jam phase. In other words, 
the phase-change between these two flows occurs with a higher density with 
the decrease of the average throughput. Furthermore, when we compare the 
three plots in Figure [H we see that with more long vehicle percentage we 
have the phase-change with a higher density. The same phenomenon is also 
shown with arrows in Figure [2j 

This phenomenon is inverted in the presence of an obstacle (on the right- 
most lan^. More precisely, when an obstacle is set we observe that the 
phase-change occurs with a lower density, see Figure [3 This phenomenon is 
more clear with the presence of long vehicles, probably because when long 
vehicles get stuck in traffic this phase emerges faster. Besides, without an 
obstacle there is an increase at the slope (absolute value) of the fundamental 
diagram where there is the wide- moving jam phase, when the percentage of 
long vehicles is increased (see Figure [2]). In other words, the flow decreases 
faster in the presence of long vehicles with the same increment of density, 
which is not observed in the case of placing an obstacle (the lines are almost 
parallel, see Figure [5|). However, this parallelism occurs only when the road 
is enough saturated. For instance in Figure [6] we see that the parallelism 
occurs in the first and second plots, but not in the third plot since the lane 
number is increased and so the emission rate per lane is decreased, i.e., there 
is not enough saturation. 

Another effect of setting an obstacle is a reduction on the traffic capacity 
as it is expected. For instance, in the experiments plotted in Figure O 
the maximum fiow that can be reached is around 1.33 veh/sec without an 
obstacle and around 1.2 veh/sec with an obstacle placed on the right-most 

^When we place an obstacle on the left-most lane, we observe that the result is almost 
the same with placing it on the right-most lane, so when we make comparisons of having 
or not having obstacle we use the result of placing it on the right-most lane. 
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Figure 4: The fundamental diagrams depending on the different average 
throughputs 




Figure 5: The fundamental and the cross-covariance diagrams without and 
with obstacle, respectively 
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Figure 6: The slope in the wide- moving jam phase with the obstacles 
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Figure 7: The effect of open road tolling on the flow phases with 3 lanes and 
4 lanes 



lane. In this figure, we also see that this reduction is more evident with the 
presence of long vehicles. 

In the case that there is an open road tolling, the wide-moving jam 
phase does not take place, so there is no anticorrelation situation. As it is 
seen in Figure [71 there is just a transition between the free flow and the 
synchronized flow. The absence of wide-moving jam phase occurs in general 
with the situations where the emission rate is low with respect to the rate 
that the vehicles are processed (the influence radius or throughput) at the 
off-toll plaza as it is expected. 

The experiments also show that the model is able to reproduce the hys- 
teresis phenomena in transition between free flow and synchronized flow 
phases (see [11]) and in transition between synchronized flow and wide- 
moving jam phases. The regions where there are saddles followed by a 
capacity drop, as in the Fig.9(e) of [6], are the regions of metastability in 
the phase-changes (from free flow to synchronized flow and from synchro- 
nized flow to wide-moving jam). This metastability phenomenon is not so 
evident when there are long vehicles. More precisely, the heterogeneity of 
traffic effects the formation of the plateaus. In Figure [2] and HI we see that in 
the case of the presence of long vehicles, the bumpy plateaus in the synchro- 
nized flow phase are replaced by more flattened plateaus, so the saddles in 
the phase-changes are not observed as they are in the absence of long vehi- 
cles. This is probably due to the fact that passenger vehicles are faster and 
so the flow of a traffic without long vehicles changes its phase more sharply. 
However, we observe that the presence of obstacles increases the number 
of the saddles even when there are long vehicles, making the fundamental 
diagrams more bumpy especially in the synchronized flow phase (see Figure 
[5] and ED. 

5 Conclusion 

It is introduced a new model for multi-lane traffic. This model is an attempt 
to define a new class of traffic models which is a hybrid between the usual 
microscopic traffic models, like car-following models, and the usual CAMs 
which we have identified as the NaSch-type models. In the process of the 
extension of the single-lane model to the multi-lane case, we have first pre- 
sented the model as an array of communicating one-dimensional CCA, and 
then we have proved that this model can be simulated by a suitable CCA. 
In this way, we have framed our multi-lane model inside the class of CCA. 
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For a first test, we have implemented the model using Pythoiij with an 
object-oriented philosophy of programming. Using a questionnaire we have 
set up two kinds of vehicles which we have used to run a series of experiments. 
Analyzing the experimental results, we have studied the influence of different 
composition of vehicles on the macroscopic behavior of the traffic in order to 
observe the typical traffic flow phenomena. We have found that the results 
are promising, indeed we reproduce and enrich the fundamental diagrams 
of traffic flow. 

Table 3: Computation Time Comparison between CPU and GPU 



Number of Vehicles 


CPU 


GPU 


(per lane) 


(sec) 


(sec) 


0; 0; 1500 


484 


45,8 


0; 0; 3000 


958 


104 


0; 0; 5000 


1608 


194 


0; 0; 10000 


4270 


556 


5000; 5000; 5000 


9679 


1288 



The code written in Python does not take advantage of CA and its 
typical synchronous behavior. For this reason, we have adapted the code 
using PyCudsl^ to partially parallelize Algorithm [2] on CPU's and we have 
seen that it is possible to boost the speed of execution by a factor of ^ 10. 
Indeed, in Table [3] we see the computation time comparison between the 
two codes made by running 1000 steps of simulation for a road with 3 lanes 
and where it is given the initial numbers of the vehicles distributed per lane 
(on a laptop equipped with a processor i7 Intel and with a graphic card 
NVIDIA CeForce CT 555M). The increased speed of execution by higher 
factors of simulation is important also to adapt the model for a forecasting 
usage. By partial parallelization, we mean that the lane-changing process 
performed using the operations ^l,^r, \ is done sequentially, and the 
only parts that can be run in parallel are the global transition functions 
'^(afe)' o.,b € {0,1}. Therefore Algorithm [2] can be completely parallelized 
only if we decide not to apply the precedence rules and assume a concurrent 
strategy of lane-changing. 

The analysis we perform here is not conclusive, but gives an insight of 
the potentiality of our model. For this reason, we suggest the following tasks 

^Thc code can bo found in the public Dropbox link 
https : //dl . dropbox . com/u/50248089/Dzsim . py , 

The code can be found in the public Dropbox link 
|https : //dl . dropbox . com/u/50248089/Cozsim . pyj 
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as future works and research directions to improve and validate the model 
and the simulator: 

• We did not use all the potentiality of the code since our aim was to 
give a first evaluation to our model. However, it would be interesting 
to consider also the experiments involving on- and off-ramps and loop- 
detectors to analyze different and more realistic situations. 

• The heterogeneity considered is reduced to two kinds of vehicles. A 
natural question is how the system reacts introducing other kinds. For 
instance, in highway environments, motorcycles or sport vehicles can 
be added to the mixed traffic of passenger vehicles and long vehicles. 
This also brings with it the interesting issue of how to tune the mem- 
bership functions for such new kinds. 

• The process of lane-changing is purely stochastic. However, in litera- 
ture there are some attempts in microscopic models where the process 
of lane-changing is described by using a fuzzy logic-based system [HH], 
thus it would be interesting to extend our model to a model in which it 
is implemented a fuzzy logic-based system to refine the lane-changing 
rules. 

• The model has to be compared with real data. In other words, a careful 
case study on specific scenarios with the data available is necessary 
for the validation by the community of people working on traffic flow 
theory, granular flow theory and traffic (transportation) engineering. 
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